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A  series  of  22  primary  reaction -control-system  engine  attitude-control  firings  were  observed  from  the  Maui  Space 
Surveillance  Site  during  the  space  shuttle  STS-115  mission.  The  firings  occurred  during  a  pass  over  Maui  on 
19  September  2006  during  which  the  orbiter  was  in  sunlight  and  the  observatory  was  in  darkness.  The  observed 
attitude  maneuvers  maintained  the  orbiter  in  an  orientation  in  which  its  long  axis  was  aligned  with  the  line  of  sight 
from  the  observatory.  This  ensured  that  the  thrust  vectors  of  all  the  observed  engine  firings  were  perpendicular  to  the 
line  of  sight,  providing  an  optimal  slde-on  observation  of  the  exhaust  The  firings  ranged  between  80  and  320  ms  in 
duration  and  involved  2  or  3  engines  for  pitch,  roll,  and  yaw  adjustments.  A  0.328  deg  field -of-vlew  acquisition  scope 
of  the  3.6  m  telescope  of  the  Advanced  Electro-Optical  System  provided  unfiltered  imagery  in  the  near-ultraviolet 
visible  spectral  region.  The  most  Interesting  white-light  features  were  transients,  one  observed  at  engine  start  up  and 
two  at  shutdown.  The  analysis  of  the  transient  speeds  reveals  that  the  startup  transient  consists  of  either  unburned 
propellant  droplets  or  higher-pressure  gas  evaporated  from  droplets  and  that  the  shutdown  transients  are 
attributable  to  a  slightly  staggered  release  of  unbumed  oxidizer  and  fuel,  respectively.  The  first  (oxidizer)  shutdown 
transient  is  the  brightest  feature,  for  which  an  intensity  evolution  analysis  is  conducted.  The  analysis  of  the  ground- 
based  data  is  fully  consistent  with  spectral  features  attributable  to  primary  reaction-control-system  engine  transients 
observed  in  previous  measurements  from  the  space  shuttle  bay  using  an  imager  spectrograph. 
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Nomenclature 

Brook  scaling  parameter  for  density,  cnr1 
Brook  scaling  parameter  for  flux,  g 
Lorentzian  scaling  parameter  for  density,  cm'1  • 
deg2 

Lorentzian  scaling  parameter  for  flux,  g  •  deg2 

Brook  angular  distribution  parameter,  unitless 

specific  heat  at  constant  pressure,  J  •  g'1  •  K'1 

particle  diameter,  pm 

initial  particle  diameter,  pm 

horizontal  video-image  frame  width  at  range,  R, 

km 

exhaust  particle  flux,  cm'2  •  s'1 
solar  spectral  flux,  W  •  cm'2 
spectral  atmospheric  transmission  and  sensitivity 
factor,  unitless 

radiant  intensity  along  image  coordinate  x 
perpendicular  to  the  thrust  axis  at  nozzle  distance 
Re,  W  *  cm-2 
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volume  emission  rate  at  nozzle  distance  Re ,  where 
y  is  the  radial  distance  with  respect  to  the  thrust 

axis,  W  cm'3 

angular  volume  emission-rate  distribution,  W  • 
cm'3 

scattered  intensity  as  a  function  of  the  solar 

scattering  angle,  W  •  sr'1 

particle  mass,  g 

number  density,  particles  cm'3 

Planck  radiation  function,  W •  sr'1  •  cm'2  •  Hz'1 

vapor  pressure,  ton- 

heat  gain  rate,  W 

Earths  hi ne  heating  rate,  W 

radiative  cooling  rate,  W 

heat  of  sublimation  cooling  rate,  W 

solar  radiation  heating  rate,  W 

range,  km 

range  vector,  km 

axial  distance  to  the  nozzle  exit,  m 
particle  radius,  pm 
temperature,  K 

orbiter  velocity  vector,  km  •  s'1 

magnitude  of  the  orbiter  velocity  vector,  km  •  s'1 

far-field  exhaust  velocity,  km  •  s'1 

projection  of  the  orbiter  velocity  vector  onto  the 

image  plane,  km  •  s'1 

image  coordinate,  m 

radial  distance  with  respect  to  the  thrust  centerline, 
m 

distance  along  the  orbiter  z  axis,  m 
elevation,  deg 

full  width  at  half-maximum  of  the  Lorentzian 
angular  distribution,  deg 
angle  between  the  orbiter  axis,  z,  and  R,  deg 
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Introduction 

THE  gaseous  and  condensed  effluents  from  spacecraft  engine 
firings  can  be  a  source  of  significant  contamination  for  sensitive 
surfaces  on  spacecraft.  This  is  particularly  problematic  in  the  case 
of  large  low- Earth- orbit  space  structures  that  require  powerful 
propulsion  systems  for  orbit  maintenance  and  attitude  control,  such 
as  the  space  shuttle  and  the  International  Space  Station  (ISS). 
Modeling  of  contaminant  fluxes  is  necessary  to  optimize  the  design 
and  integration  of  spacecraft  subsystems  and  instruments.  The 
models  require  precise  information  on  exhaust  constituents  and  their 
angular  distribution.  Condensed  exhaust  matter  is  a  primary  concern 
because  of  the  high  velocity  of  exhaust  impacts  and  the  chemically 
sensitive  nature  of  optical  surfaces 

Instruments  on  the  ISS  are  exposed  to  contamination  by  service 
vehicles,  including  the  space  shuttle,  Progress,  Soyuz,  and  the 
ESA  automated  transfer  vehicle  during  docking  and  separation 
maneuvers  A  series  of  experiments  have  investigated  contamination 
caused  by  space  shuttle  reaction-control -system  engine  exhaust 
[1,2].  The  experiments  measured  deposition  rates  and  clearly 
identified  significant  surface  erosion,  as  evidenced  by  microcrater 
formation  attributed  to  hypervelocity  impacts  by  condensates  in  the 
exhaust.  The  measurements  were  only  conducted  on  the  thrust 
centerline,  and  therefore  the  deposition  rates  from  these  bipropellant 
engines  as  a  function  of  angle  with  respect  to  the  thrust  axis  are  not 
known. 

The  STS- 115  mission  in  September  2006  provided  a  valuable 
opportunity  to  examine  condensed  exhaust  species,  including 
angular  distributions,  during  a  terminator  pass  over  the  Maui  Space 
Surveillance  Site  (MSSS).  A  set  of  dedicated  bums  was  planned  as 
part  of  the  U.S.  Air  Force  Research  Laboratory  Maui  Analysis  of 
Upper- Atmospheric  Injections  (MAUI)  experiment.  Because  the 
unique  lighting  conditions  (orbiter  in  sunlight  and  observatory  in 
darloiess)  were  optimal  for  imaging,  and  because  there  was  concern 
over  the  possibility  of  having  lost  the  parachute  hatch,  the  flight 
officer  canceled  the  dedicated  MAUI  bums  and  ordered  an  imaging 
overflight  during  which  the  shuttle’s  long  axis  was  aligned  with  the 
line  of  sight  of  the  MSSS  telescope  (landmark  mode).  Although 
unintended,  the  requirement  to  maintain  a  fixed  attitude  with  respect 
to  the  observatory  resulted  in  a  long  series  of  attitude-control  firing 
bursts.  This  provided  an  unprecedented  opportunity  to  acquire  video 
imagery  of  sunlight  scattering  by  exhaust  condensed  matter  during  a 
large  number  of  engine  startups  and  shutdowns. 

The  purpose  of  the  present  study  was  to  denve  condensate  speeds 
and  angular  fluxes  from  remote  observations  of  sunlight  scattering 
from  transient  condensate  clouds.  The  current  paper  presents  an 
analysis  of  transients  observed  for  a  total  of  22  attitude-control  bums 
involving  2  or  3  primary  reaction-control-system  (PRCS)  engines. 
These  are  the  most  frequently  used  space  shuttle  engines  during  ISS 
proximity.  For  each  bum,  3  transients  were  identified,  and  the  speeds 
of  the  center  of  the  observed  features  could  be  accurately  determined 
The  retrieved  speeds  allowed  the  determination  of  the  ejection  times 
of  the  matter  causing  the  transients  and  their  association  with  engine 
startup  or  shutdown.  For  the  first  time,  angular  distributions 
determined  from  the  transverse  intensity  profiles  of  the  transient 
clouds  are  reported.  We  have  also  conducted  an  analysis  of  the 
evolution  of  the  integrated  intensity  of  the  most  intense  transient 
for  each  bum,  which  we  attribute  to  pure  oxidizer  condensates. 
The  observations  are  compared  with  model  calculations  of 
Mie- scattering  intensities  from  condensed  sublimating  propellant 


particles.  We  also  relate  the  present  observations  to  previously 
unreported  spectra  associated  with  particle  scattering  observed 
following  engine  bums  in  sunlight  using  the  Arizona  Airglow 
Instrument  (GLO)  [3,4]  imager  spectrograph  flown  in  the  shuttle  bay 
during  STS-63. 


Experimental 

STS-l  15  on  Atlantis  was  launched  on  9  September  2006.  The 
observation  pass  occurred  on  19  September  2006  following  un¬ 
docking  from  the  ISS,  Culmination  with  respect  to  the  observatory  at 
MSSS  was  at  15:29:48  h  coordinated  universal  time  (UTC).  Figure  1 
provides  a  summary  of  pass  parameters  as  a  function  of  time, 
including  range  R,  elevation  or,  solar  phase  angle  <p ,  and  the  angle  y 
between  the  orbiter  axis  z  and  R.  The  geometry  of  the  experiment 
and  associated  parameters  are  shown  schematically  in  Fig.  2.  The 
space  shuttle  flew  over  Hawaii  in  a  landmark  mode  in  which  the 
orbiter  z  axis  was  oriented  parallel  to  the  line  of  sight  such  that  high- 
resolution  imagery  could  be  acquired  of  the  shuttle  aft  section.  The 
bums  that  are  the  subject  of  this  work  were  a  series  of  attitude 
maneuvers  fired  perpendicularly  to  the  orbiter  z  axis  to  maintain 
landmark  mode.  The  bottom  chart  in  Fig,  1  provides  markers  at  the 
time  of  those  attitude-control  maneuvers.  The  time  interval  during 


Fig.  1  Observational  parameters  during  the  19  September  2007  pass  of 
Atlantis  over  Maui. 
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Fig.  2  Schematic  representation  of  the  observation  geometry. 


which  quality  imaging  of  bums  could  be  acquired  is  indicated  with 
vertical  dashed  lines  After  15.29:40  h,  cirrus  clouds  caused  signif¬ 
icant  background  radiation  as  MSSS  approached  sunrise,  substan¬ 
tially  reducing  the  quality  of  the  observations.  The  orbiter  velocities 
in  the  Earth-centered-inertial  (ECI)  and  Earth-centered- fixed  co¬ 
ordinate  systems  were  7.715  and  7.42  km  •  s_1,  respectively  The 
orbital  altitude  was  335.15  ±  0.20  km  within  the  time  interval  of 
the  investigated  bums,  with  R  declining  from  780  to  446  km  dunng 
the  approach.  As  shown  in  Fig.  1,  during  the  investigated  bum 
interval,  the  angle  /  remained  nearly  constant  at  a  value  of  176.3  deg 
to  within  0.5  deg.  The  minor  discrepancy  with  a  perfect  alignment 
has  minimal  bearing  on  the  present  analysis. 

Wide-field-of-view  (FOV)  imagery  of  the  space  shuttle  main¬ 
taining  axial  alignment  with  the  MSSS  line  of  sight  was  obtained 
using  the  large-aperture  acquisition  telescope  (LAAT).  This  0.58-m- 
aperture  telescope,  situated  on  the  advanced  electro-optical  system 
mount,  was  operated  with  a  nominal  FOV  of  0.328  deg.  The 
lntensified-charge-coupled-device  (1CCD)  detector  exhibits  sensi¬ 
tivity  in  the  spectral  range  of  0.4  to  0  9  fim.  The  sensor  gain  was 
set  to  limit  the  size  of  the  spot  attributable  to  saturation  by  the 
solar-illuminated  orbiter  hard  body.  Unfortunately,  a  sensitivity 
calibration  was  only  conducted  for  a  configuration  involving  filtered 
imagery  with  a  different  gain  setting.  The  change  to  imaging 
experiments  came  so  late  that  calibration  measurements  could  not 
be  conducted  for  the  reduced-gain  setting  of  the  LAAT  sensor. 

Table  1  lists  the  pertinent  parameters  for  22  attitude -control  bums 
for  which  quality  (30  frames  per  second)  video  footage  was  obtained. 
Most  of  the  bums  were  the  minimum  bum  length  of  0.08  s.  Bums  2 
and  8  were  0. 16  s  long,  and  bum  9  was  0.24  s  long.  Three  types  of 
attitude -control  maneuvers  actuating  2  or  3  PRCS  engines  were 


conducted  pitch  down,  roll  left,  and  yaw  right  The  F1U,  L3D,  and 
R3D  engines  were  used  for  pitch  down;  the  LIU  and  R3D  engines 
were  used  for  roll  left;  and  the  F3L  and  R1R  engines  were  used  for 
yaw  right.  The  thrust  axis  for  all  of  these  engines  is  perpendicular  to 
the  shuttle  axis  and  is  therefore  essentially  perpendicular  to  the  line  of 
sight,  providing  a  side -on  view  of  the  engine  exhaust.  A  space  shuttle 
PRCS  engine  has  a  nominal  thrust  of  3880  N  (870  lbf)  and  uses  a 
monomethyl-hydrazine  (MMH)  fuel  and  N204  oxidizer. 

Table  1  also  lists  the  projection  of  the  orbiter  ECI  velocity 
vector  v  onto  the  plane  normal  to  the  line  of  sight.  This  projection  is 
directly  proportional  to  the  angular  velocity  of  stars  streaking  across 
the  FOV  as  the  orbiter  is  tracked  and  is  given  by 

v±  =  vsin(Z{R,  v})  (1) 

The  dimensions  of  the  video  frame  at  the  orbiter  altitude  can 
therefore  be  calibrated  by  determining  the  velocity  of  the  star  streaks 
in  units  of  pixels  per  unit  time,  from  which  the  scaling  factor 
(distance  per  pixel)  can  be  derived.  This  approach  neglects  a 
relatively  insignificant  contribution  due  to  the  motion  of  the  Earth 
At  the  small  nominal  FOV  of  0.328  deg,  it  is  within  the  accuracy  of 
the  present  analysis  to  assume  tan(FOV)  ^  FOV,  thereby  allowing 
the  use  of  planar  geometry  for  the  image  analysis.  Table  1  lists  the 
derived  horizontal  frame  widths  d ,  which  vary  from  ~3-5  km, 
versus  orbiter  distance.  From  FOV  =  d/R ,  we  derive  an  average 
FOV  of  0.374  deg  for  the  22  bums  with  a  small  standard  deviation 
of  0.004  deg,  thereby  confirming  a  consistent  calibration  of  the 
observed  video  imagery. 

The  present  analysis  of  STS- 1 1 5  observations  is  complemented  by 
spectral  measurements  conducted  with  the  GLO  imager  spectro¬ 
graph  flown  earlier  on  STS-63.  The  GLO  spectral  measurements  are 
reported  here  for  the  first  time.  The  spacebome  GLO  instrument  has 
been  described  in  detail  before  [3,4],  including  GLO  measurements 
of  engine-exhaust  interactions  with  the  upper  atmosphere  reported 
by  Bernstein  et  al.  [5,6].  GLO  was  mounted  in  the  rear  section  of 
the  shuttle  bay.  It  comprises  a  nine-section  spectrograph,  three 
monochromatic  imagers,  and  a  TV  camera,  all  boresighted  to  view  in 
the  same  direction.  The  spectrograph  and  imagers  have  1CCD  focal  - 
plane  detectors.  The  nine  slightly  overlapping  spectrograph  sections 
permit  simultaneous  recording  of  the  spectrum  from  115.0  to 
900.0  nm  with  a  spectral  resolution  of  about  4  to  10  A.  Their 
combined  focal-plane  image  is  4500  pixels  wide  in  the  wavelength 
(dispersion)  dimension,  perpendicular  to  the  slit,  and  192  pixels  in 
the  spatial  dimension,  along  the  slit.  The  slit  image  at  the  detector  is  a 
narrow  portion  of  the  image  of  the  distant  object  being  observed, 
preserving  spatial  resolution  in  the  slit-length  direction  The  design 


Table  1  List  of  PRCS  engine  burns  analyzed  in  this  work 


Bum  no 

UTC 

Maneuver 

Bum  length,  s 

Range,  km 

Ujl  ,  km/s 

d,  km 

1 

15:28:18  993 

Pitch  down 

0.08 

778.837 

4.607 

4.976 

2 

15:28  21,073 

Yaw  right 

0  16 

766.504 

4.670 

4.924 

3 

15:28.25.473 

Pitch  down 

008 

740.743 

4  809 

4.770 

4 

15:28:44.593 

Pitch  down 

0.08 

635.166 

5.508 

4.137 

5 

15:28  45.313 

Pitch  down 

0.08 

631,436 

5.537 

4.113 

6 

15:28  48.353 

Yaw  right 

008 

615  922 

5.663 

3,990 

7 

15:28:51.073 

Yaw  right 

0.12 

602.383 

5.779 

3.914 

8 

15.2900.673 

Pitch  down 

0.16 

557.521 

6.203 

3.665 

9 

15:2902.673 

Roll  left 

024 

548.819 

6.294 

3.567 

10 

15:29:03.073 

Roll  left 

0.08 

547.108 

6.312 

3.551 

11 

15:29:08.673 

Pitch  down 

0.08 

524.242 

6.568 

3.442 

12 

15:29:12.193 

Pitch  down 

0.08 

510.990 

6.726 

3.342 

13 

15:29:14.113 

Roll  left 

0.08 

504.157 

6.812 

3.257 

14 

15:29:16.033 

Roll  left 

0.08 

497.619 

6.896 

3.251 

15 

15:29:16.593 

Yaw  nght 

0.08 

495.768 

6920 

3.250 

16 

15:29:18.993 

Pitch  down 

0.08 

488.141 

7.021 

3.181 

17 

15.29:23.713 

Pitch  down 

0.08 

474.638 

7.209 

3.105 

18 

15:29:26.193 

Roll  left 

0.08 

468.385 

7  300 

3.067 

19 

15:29:27.553 

Roll  left 

0.08 

465.215 

7.347 

3.088 

20 

15:29:29.473 

Pitch  down 

0.08 

461.060 

7.410 

3  047 

21 

15:29:38.313 

Pitch  down 

0.08 

447.029 

7.630 

2.940 

22 

1529:39.073 

Roll  left 

0.08 

446.231 

7.643 

2.833 
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of  the  slit  and  fore  optics  affords  a  field  of  view  of  0.2  x  8.5  deg.  In 
the  GLO  measurements,  the  instrument  was  pointed  almost  parallel 
to  the  thrust  axis  of  a  PRCS  engine  This  orientation  allowed  the 
observation  of  far-field  exhaust  radiance  while  avoiding  radiance  due 
to  the  hot  near-field  vacuum  core.  Data  acquisition  and  downloading 
are  described  in  separate  publications  [7,8]. 

Results 

Close  inspection  of  the  video  imagery  leads  to  the  observation 
that  3  transients  accompany  every  bum;  this  is  demonstrated  for 
bum  17  in  Fig.  3.  As  we  will  show,  the  first  transient  is  associated 
with  engine  startup,  and  the  second  and  third  transients  are  associated 
with  engine  shutdown.  At  the  time  of  bum  17,  the  horizontal  width  of 
the  frame  corresponded  to  3. 1 05  km  at  the  range  of  the  orbiter,  which 
is  at  the  center  of  the  circular  region  of  saturated  pixels.  Stray  light 
due  to  the  intense  sunlight  scattered  off  the  body  of  Atlantis  caused 
significant  blooming.  The  sharp  streaks  are  stars,  one  of  which  was 
used  for  the  frame-width  calibration.  Bum  17  was  a  pitch-down  bum, 
which  involves  the  two  aft  engines  L3D  and  R3D  firing  downward 
and  the  forward  engine  F1U  firing  upward  (toward  the  bottom  right 
comer  of  the  image)  in  the  shuttle  frame  of  reference.  The  first 
transient  is  the  dimmest  and  has  the  highest  velocity  with  respect  to 
the  shuttle.  The  second  transient  is  the  brightest  and  most  compact, 
and  the  third  has  the  lowest  velocity  and  forms  a  trail  due  to  a  slow 
tail-off  in  the  expelled  material  following  the  initial  burst.  The 
imagery  also  shows  that  there  is  considerable  interaction  between  the 
visible  exhaust  features  of  the  two  aft  engines 
The  velocities  of  the  individual  transients  were  analyzed  for  every 
bum  and  engine  by  recording  the  time  dependence  of  the  distance 
of  the  point  of  highest  intensity  of  the  respective  feature  from  the 
center  of  the  orbiter  along  the  exhaust  centerline.  A  frame-by-frame 
analysis  established  that  the  transient  velocities  are  constant  with 
respect  to  the  distance  of  the  visual  center  of  the  feature  to  the  orbiter, 
Thus,  there  is  no  effect  due  to  interaction  with  the  upper  atmosphere, 
which  is  not  surprising  given  the  low  thermospheric  densities  at 
330  km  during  the  near  solar-minimum  conditions  of  the  experiment 
The  analysis  of  the  velocities  is  summarized  in  Table  2,  in  which 
the  average  transient  velocities  are  sorted  by  engine  and  transient. 
The  global  averages  of  transient  velocities  averaged  over  all  22  bums 
are  also  provided.  There  are  no  noteworthy  differences  between  the 
engines,  the  respective  average  values  lying  close  to  the  standard 
deviations  of  the  global  average  for  a  respective  transient.  The  first 
transient  has  an  average  velocity  of  3.13  ±  0 .19  km  •  s_1,  close  to 
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Fig.  3  Video  frames  separated  by  0.2  s  of  bum  17. 


Table  2  Axial  transient  velocities  in  km  •  s  1  determined  for  bums 
listed  in  Table  1. 


Engine 

Transient  1 

Transient  2 

Transient  3 

F1U 

3.268 

± 

0.105 

1  599 

± 

0.067 

0.795 

± 

0.093 

L3D 

3.190 

± 

0  172 

1.686 

± 

0.225 

0.655 

± 

0.143 

R3D 

3.078 

± 

0  170 

1.679 

± 

0.172 

0.701 

± 

0.152 

F3L 

3.093 

± 

0.175 

1.485 

± 

0,097 

0.631 

± 

0.051 

R1R 

3.348 

± 

0.080 

1.392 

± 

0.146 

0.675 

± 

0.090 

LIU 

2.941 

± 

0.205 

1.527 

± 

0.086 

0.713 

± 

0.033 

Transient  Average 

3.132 

± 

0.185 

1.630 

± 

0.175 

0.707 

± 

0.128 

the  terminal  exhaust  velocity  of  a  PRCS  engine  [9,10].  The  average 
velocity  of  transient  2  of  1.63  ±  0.18  km  •  s-1  is  approximately  half 
the  velocity  of  transient  1,  and  transient  3  has  an  average  velocity  of 
0.71  dr  0  13  km  s-1,  approximately  half  that  of  transient  2.  No 
noticeable  trend  with  time  of  bum  is  observed,  and  therefore  effects 
due  to  engine  heating  on  the  terminal  velocities  can  be  excluded 
Having  derived  the  average  velocities  of  the  transients,  the  precise 
time  of  origin  of  the  transients  could  be  determined  The  times  of 
origin  of  transient  1  are  identical  to  the  engine  startup  time,  and  the 
times  of  origin  of  transients  2  and  3  are  essentially  the  same  and  occur 
later  The  intervals  between  the  origin  times  of  transient  2  and  3  and 
engine  startup  agree  closely  with  the  bum  times  reported  to  us  by 
NASA.  Consequently,  transients  2  and  3  can  be  finked  to  engine 
shutdown 

Figure  4  shows  intensity  profiles  through  the  transient  images  of 
the  F1U  engine  shown  in  Fig.  3  Each  pixel  unit  corresponds  to 
6.47  m.  Profiles  transverse  to  the  thrust  axis  provide  information  on 
the  angular  flow  distribution  of  the  respective  transient.  The  top 
figure  is  a  transverse  profile  of  transient  2  at  the  time  of  the  center 
frame  in  Fig.  3.  The  data  are  compared  with  a  Gaussian  fit  to  the  data, 
allowing  for  an  offset  with  respect  to  background  and  maximum.  It  is 
seen  that  a  Gaussian  is  a  good  representation  of  the  distribution.  In 
the  central  plot  of  Fig.  4,  a  transverse  profile  of  transient  3  is  shown, 
as  observed  in  the  bottom  frame  of  Fig.  3,  In  this  case,  a  Lorentzian 
provides  a  significantly  improved  fit.  The  bottom  plot  in  Fig.  4  is  a 
longitudinal  profile  taken  from  the  bottom  frame  in  Fig.  3,  running 
from  the  orbiter  position  along  the  F1U  thrust  centerline.  The 
longitudinal  plot  has  been  corrected  with  respect  to  background,  as 
derived  from  the  fits  in  the  top  2  plots  of  Fig  4.  The  background 
is  consistently  found  to  be  approximately  50  intensity  units.  The 
intensity  scale  spans  256  units,  and  a  nonlinear  intensity  regime  must 
be  assumed  as  the  maximum  intensity  is  approached  The  orbiter 
presents  the  most  intense  saturated  feature  at  around  200  intensity 
units.  Both  the  second  and  third  transients  of  the  FI  U  engine  can  be 
identified  as  peaks  in  the  profile.  It  is  also  apparent  that  there  is  some 
overlap  between  the  two  features. 

The  angular  distribution  of  the  exhaust  transient  constituents  are 
obtained  from  the  images  in  Fig.  3  by  applying  an  Abel  transfor¬ 
mation  [  1 1  ]  to  a  series  of  slices  through  the  transient  cloud  transverse 
to  the  thrust  axis  at  equally  spaced  axial  distances  to  the  nozzle 
exit,  R(.  We  assume  that  the  feature  is  optically  thin  and  axially 
symmetric.  The  transformation  converts  the  relative  radiant  line-of- 
sight  intensity  distributions  /(/?,,  x),  where  x  is  the  image  coordinate 
transverse  to  the  thrust  axis,  to  radial  volume  emission-rate 
distributions  /(/?,,  y),  where  y  is  the  radial  distance  with  respect  to 
the  thrust  centerline  The  angular  volume  emission-rate  distribution 
1(0)  averaged  over  the  Re  interval  of  the  selected  slices  is  then 
obtained  by  using  an  angular  binning  approach  across  all  points 
/(/?,,  y)  of  the  Abel-transformed  cloud  with  the  relationship 

6=  tan-'ty/fl,)  (2) 

Figure  5  shows  the  derived  angular  distributions  for  the  first 
and  second  transients  In  the  top  chart,  the  second  transient  angular 
distribution  is  compared  with  the  exhaust-gas  angular  density 
distribution  derived  from  the  expression  for  a  Brook  source  [12], 


N(Re,  0)  =  A  exp[— £(1  -  cos(0))]/?72 


(3) 
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Fig.  4  Transverse  and  longitudinal  intensity  profiles  of  bum  17 
transients  during  indicated  times. 

using  B  =  8.06  specified  for  a  PRCS  engine.  The  density  is  pro¬ 
portional  to  the  emission  rate  at  optically  thin  conditions.  The 
comparison  indicates  that  the  derived  angular  distribution  is  slightly 
narrower  than  that  specified  for  the  PRCS  exhaust  gases.  Also  shown 
is  a  fit  of  the  Brook  parameter  B  to  the  observed  distribution,  yielding 
B  =  14  4  for  the  transient. 

The  bottom  chart  in  Fig.  5  compares  the  derived  angular  distri¬ 
bution  of  transient  3  to  3  calculated  distributions:  1 )  the  PRCS  Brook 
distribution,  2)  an  attempt  to  fit  a  Brook  distribution  to  the  data,  and 
3)  a  Lorentzian  fit.  The  Lorentzian  is  given  by  the  expression 

=  +  (4) 

The  angular  parameters  B  and  T  used  for  the  calculated  distributions 
in  Fig.  5  are  listed  in  Table  3,  along  with  the  Brook  values  determined 
for  transient  2  and  those  applicable  to  the  gaseous  exhaust  of  a  PRCS 
engine.  Although  a  Brook  distribution  provides  an  adequate 
description  at  angles  below  30  deg,  at  larger  angles,  there  is  clearly 
additional  intensity,  suggesting  that  the  angular  distribution  of 
transient  3  does  not  obey  a  Brook  distribution.  The  Lorentzian 
reproduces  the  wide-angle  signal  well  in  the  I(x)  distributions 
(middle  frame  of  Fig.  4);  however,  a  fit  of  a  Lorentzian  to  the 
transformed  1(0)  does  not  capture  the  distribution  correctly  at  the 
highest  angles  (bottom  frame  of  Fig.  5). 

To  derive  more  information  on  the  source  of  the  transient 
signatures,  the  temporal  evolution  of  the  spatially  integrated  intensity 
of  an  entire  transient  is  examined.  This  was  conducted  for  transient  2 
because  it  is  the  most  intense  and  compact.  Figure  6  shows  the 


Angle  (deg) 

Fig.  5  Comparison  of  derived  transient  angular  distributions  with 
PRCS  Brook  flow  and  parametric  fits. 


dependence  of  the  integrated  intensity  of  the  second  transient  as  a 
function  of  time  following  engine  shutdown  for  the  F1U  engine 
in  pitch-down  maneuvers  1  and  17  The  intensity  was  determined 
following  background  subtraction  by  defining  the  image  area  of 
the  transient  with  a  polygon.  Critical  in  this  effort  is  the  precise 
determination  of  background  intensities,  because  errors  are  additive 
in  the  intensity  integration.  The  intensities  are  normalized  with 
respect  to  the  range  (R~2)  and  atmospheric  transmission  In  the  case 
of  the  latter,  an  average  wavelength  of  600  nm  was  assumed  After 
~0.5  s,  the  integrated  intensity  declines  with  time.  Before  0.5  s,  the 
intensity  appears  to  be  constant  or  rising  with  time.  Examination  of 
the  individual  pixel  intensities  indicates  that  in  the  case  of  bum  1 ,  no 
pixels  registered  intensities  exceeding  200  units  before  background 
subtraction,  whereas  in  the  case  of  bum  17,  no  pixels  had  intensities 
exceeding  210  units  beyond  0.3  s.  It  can  therefore  be  concluded 
that  saturation  and  nonlinearity  in  the  charge-coupled-device  (CCD) 
response  are  most  likely  not  the  source  of  the  varying  time  depen¬ 
dence  of  the  integrated  intensity  of  the  second  transient 
The  GLO-2  mission  on  STS-63  that  flew  in  February  1995 
obtained  spectral  data  on  PRCS  engine  firings  during  daytime 
conditions  at  similarly  low  thermospheric  densities  at  an  orbital 
altitude  of  nearly  400  km  [5].  Figure  7  shows  spectra  recorded 
during  a  3  s  PRCS  bum  and  approximately  3  s  after  the  bum.  The 
spectrograph  line  of  sight  was  nearly  parallel  to  the  thrust  centerline, 
and  the  spectral  acquisition  times  were  2  s.  The  spectral  features  in 
the  bottom  spectrum,  which  is  on  an  expanded  vertical  scale  in 
comparison  with  the  top  spectrum,  are  due  to  OH,  NO,  and  NH 


Table  3  Density  (cm-3)  and  deposition  (g  •  cm-2)  parameters  derived 
from  the  angular  distribution  analysis  of  the  shutdown  transients 


Density 

Deposition 

A,  cm-1  B 

r,  deg 

*.g 

A\  g  •  deg2 

Gaseous  exhaust 

1.31  x  1020  8  07 

_ 

_ 

_ 

Transient  2 

-  14.4 

— 

1.12 

— 

Transient  3 

8.95 

35 

1.12 

340 
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Fig.  6  Range  and  transmission  normalized  intensity  versus  time  after 
shutdown  for  2  F1U  second  transients. 


Wavelength  (nm) 

Fig.  7  GLO  spectra  of  exhaust  far-field  radiance  during  (/  =  0)  and 
after  a  PRCS  bum  in  daylight. 

exhaust  species  (see  [5]  for  an  assignment  of  the  bands).  The  top 
spectrum  exhibits  abroad  background  with  increasing  intensity  with 
wavelength.  The  broad  feature  extends  to  longer  wavelengths.  The 
two  sharp  features  at  —393  and  397  nm  are  the  well-known  Ca  II  K 
and  H  lines  observed  in  the  solar  spectrum.  The  intense,  otherwise 
featureless,  band  can  therefore  be  primarily  attributed  to  solar 
scattering. 

Figure  8  shows  the  time  evolution  of  the  integrated  intensity  of  the 
GLO  data  between  430  and  460  nm  for  the  same  bum  that  produced 
the  data  in  Fig.  7,  identified  as  PRCS  ram,  and  an  additional  PRCS 
bum,  identified  as  PRCS90.  The  ram  and  90  designations  refer  to  the 
angle  of  attack,  ram  being  a  bum  in  which  the  exhaust  vector  points 
into  the  direction  of  the  velocity  vector,  and  90  refers  to  a  cross  bum. 
The  radiance  in  Fig.  8  is  on  an  absolute  scale, 

\R  =  \  rayleigh  =  106  photons  cm-2  •  s'  1  radiated  into  sr 

and  the  time  corresponds  to  an  estimate  of  the  end  of  the  signal 
acquisition  interval.  In  the  case  of  the  cross  bum  (PRCS90),  the  first 
data  point  is  associated  with  a  spectrum  that  reveals  both  exhaust 


Fig.  8  Comparison  between  measured  (GLO)  and  calculated  (1  /im 
oxidizer  droplet)  430-460  nm  band  intensity  evolution. 


spectral  features  and  the  significantly  more  intense  continuous 
radiation  The  interval  therefore  overlaps  with  t  —  0 

The  two  bums  at  different  angles  of  attack  exhibit  very  similar 
radiances  after  shutdown.  Given  the  high  altitude,  this  is  expected,  as 
effects  due  to  atmospheric  interactions,  which  are  sensitive  to  the 
thrust  angle  of  attack,  are  minimized.  The  intensity  rapidly  dims  in 
the  first  seconds  following  the  end  of  the  bum  and  then  declines  at  a 
significantly  lower  rate.  This  observation  is  consistent  with  the 
interpretation  that  solar  scattering  is  the  source  of  transients  2  and  3 
that  we  associate  with  engine  shutdown.  The  sharp  drop  in  intensity 
can  be  related  to  transient  2,  which  recedes  from  the  orbiter  with  a 
higher  velocity,  and  the  fact  that  the  scattering  intensity  recorded 
with  the  GLO  instrument  scales  with  R~2.  The  slow  decline  must  be 
associated  with  the  persistent  trails  of  the  slower  transient  3.  The 
solid  line  in  Fig.  8  is  the  calculated  intensity  evolution  of  a  narrow 
pulse  of  l-/xm-diam  droplets  flying  from  the  GLO  sensor  at  the 
transient  2  representative  speed  of  1 .6  km  •  s_1 .  The  model  assumes 
that  transient  2  is  due  to  solar  scattering  from  oxidizer  condensates,  as 
will  be  discussed  in  the  next  section.  The  model  results  are  consistent 
with  the  early  fast  decline. 

Discussion 

Qualitative  Interpretation  of  Transients 

The  combined  ground-  and  space-based  observations  of  PRCS 
engine  transients  in  sunlight  suggest  that  the  transients  arise  from 
solar  scattering  from  condensates  in  the  exhaust.  The  present  data 
therefore  provide  an  alternative  approach  to  deposition  measure¬ 
ments  in  examining  contamination  associated  with  shuttle  PRCS 
engine  firings.  Solar  particle  scattering  in  association  with  space 
shuttle  experiments  has  been  previously  observed  in  a  water  vent 
experiment  that  occurred  during  a  similar  terminator  pass  over 
MSSS  [13].  In  that  case,  the  far-field  solar  scatter  was  primarily 
associated  with  scattering  from  submicron -sized  particles  formed  by 
recondensation  of  water  sublimated  from  larger  ice  particles  emitted 
at  the  orifice.  The  high  velocity  of  transient  1  formed  at  engine  startup 
suggests  that  this  transient  may  be  attributable  to  condensation  of 
some  exhaust  species. 

Condensation  of  exhaust  is  typically  observed  for  liquid-oxygen/ 
liquid-hydrogen  boosters,  as  first  reported  by  Kung  et  al.  [14]. 
Exhaust  condensation,  however,  would  imply  that  the  solar  scatter 
would  be  observable  throughout  the  bum,  which  is  not  the  case,  as 
demonstrated  by  the  low  background  levels  in  the  bottom  spectrum 
of  Fig.  7  recorded  at  steady-state  bum  conditions.  Alternatively,  it  is 
conceivable  that  at  engine  startup,  before  reaching  steady  state, 
mixing  of  oxidizer  and  fuel  are  not  complete,  leading  to  unbumed 
droplets  that  take  on  the  velocity  of  the  exhaust  through  collisional 
drag  forces.  Thus,  it  is  plausible  that  condensates  reside  at  the  exhaust 
front.  Droplets  at  the  exhaust  front  can  also  lead  to  a  higher-gas- 
density  region  due  to  evaporation  of  the  droplets  in  flight.  The  higher- 
density  gas  is  then  observed  through  enhanced  solar-induced 
fluorescence  intensities.  Solar-induced  fluorescence  should  have 
high  yields  in  the  case  of  N02,  the  primary  evaporation  product  of 
N204,  given  the  significant  absorption  cross  sections  of  N02  in  the 
visible  spectrum  Meanwhile,  the  lower  sublimation  rates  and  higher 
heats  of  sublimation  of  MMH  point  to  a  higher  probability  of 
condensate  survival  for  the  fuel  versus  the  oxidizer 

It  is  worth  noting  that  the  first  transient  appears  more  compact  and 
brighter  in  the  case  of  the  LIU  engine  used  during  the  roll-left 
maneuvers  The  thrust  axis  of  this  engine  on  the  aft  left  pod  points  in 
the  direction  of  the  vertical  stabilizer  of  the  orbiter.  It  is  possible  that 
interaction  between  the  exhaust  and  the  stabilizer  surfaces  are  the 
cause  of  increased  brightness  of  the  first  transient  of  this  engine.  For 
all  other  engines  used  in  the  present  experiments,  the  exhaust  does 
not  interact  directly  with  orbiter  surfaces. 

The  shutdown  transients  are  most  likely  due  to  small  amounts  of 
oxidizer  and  fuel  present  in  the  vacuum  manifold  between  the  tank 
valves  and  the  combustion  chamber.  The  manifold  volume  of  a  space 
shuttle  PRCS  engine  is  cm-3  for  the  oxidizer  and  ~8  cm”3 
for  the  fuel.  When  the  oxidizer  and  fuel  valves  are  turned  off 
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simultaneously,  the  oxidizer  and  fuel  in  the  vacuum  manifold  are  no 
longer  pressurized.  They  escape  the  manifold  through  their  vapor 
pressure  The  more  volatile  oxidizer  can  then  be  expected  to  vent 
into  the  combustion  chamber  at  higher  rates  than  the  fuel,  which  is 
also  considerably  more  viscous.  This  interpretation  implies  that 
transient  2  is  due  to  oxidizer  condensates  and  transient  3  is 
attributable  to  the  MMH  fuel. 

The  size  of  the  condensates  dictates  both  the  scattering  intensities 
and  transient  velocities.  If  all  transients  are  due  to  condensates,  their 
terminal  velocities,  listed  in  Table  2,  are  governed  by  their  interaction 
with  the  expanding  gaseous  species.  At  a  low  size  limit,  particles  are 
expected  to  take  on  the  terminal  velocity  of  the  gas.  This  is 
demonstrated  in  Fig.  9,  which  plots  the  particle  velocity  as  a  function 
of  distance  from  the  nozzle  for  3  particle  sizes.  The  calculations 
assume  that  particles  are  expanded  with  the  steady-state  exhaust 
of  a  PRCS  engine  at  the  exit  (—3  km  s_l)  and  that  full  gas 
accommodation  occurs  (i.e.,  every  gas-particle  collision  leads  to 
sticking  of  the  gas  species).  This  corresponds  to  the  conditions  for 
maximum  acceleration,  and  thus  the  determined  terminal  velocities 
must  be  regarded  as  upper  limits.  The  model  determines  that 
condensates  of  a  0  1  /on  diameter  experience  no  velocity  slip, 
whereas  1-  and  10-^m-diam  particles  reach  asymptotic  velocities  of 
2850  and  2100  m/s,  respectively.  In  the  case  of  the  startup  transient 
(first  transient),  the  condensates  would  have  essentially  the  terminal 
velocity  of  the  exhaust  gases  in  the  collisionless  regime.  Thus,  the 
calculations  in  Fig.  9  suggest  that  if  the  first  transient  is  due  to 
condensates,  the  average  particle/droplet  size  is  significantly  below 
1  /*m.  A  subsequent  discussion  of  the  Mie-scattering  properties 
demonstrates  that  the  scattering  cross  sections  decline  rapidly  with 
particle  size  in  this  size  regime,  which  would  be  consistent  with 
the  low  intensity.  However,  given  the  available  information,  it 
cannot  be  definitively  concluded  whether  the  first  transient  is  due  to 
fluorescence  from  a  higher-gas-density  region  or  direct  scattering 
from  condensate  particles.  The  calculations  in  Fig.  9  do  not  apply 
directly  to  the  shutdown  transients  (transients  2  and  3)  because  the 
condensates  are  not  entrained  in  the  same  steady-state  PRCS 
exhaust-gas  flow. 

We  can  interpret  the  second  transient  as  being  due  to  submicron- 
particle  condensates  In  this  interpretation,  a  very  lean  propellant 
mixture  at  shutdown  leads  to  a  lower-temperature  expansion  and 
terminal  velocities  of  —  1.6  km*  s_1.  The  exhaust  gases  resulting 
from  the  leaner  combustion  would  be  expected  to  accelerate  the 
oxidizer  droplets  to  lower  velocities  than  those  produced  at  optimal 
combustion  conditions.  Alternatively,  the  velocities  can  be  ex¬ 
plained  by  complete  and  instantaneous  vaporization  of  the  liquid 
manifold  oxidizer  when  striking  the  hot  combustion-chamber  walls. 
Given  the  —7  6  g  mass  of  the  manifold  oxidizer  and  assuming  a 
density  of  1.52  g  •  cm- 3  and  a  chamber  volume  of  —320  cm-3, 
complete  vaporization  should  produce  chamber  pressures  of  several 


Fig.  9  Particle  velocity  as  a  function  of  axial  nozzle  distance  for  3 
particle  sizes. 


atmospheres.  Thus,  the  transient  may  be  associated  with  a  hard 
expansion  of  purely  gaseous  oxidizer  that  leads  to  recondensation. 
In  such  an  expansion,  the  thermal  energy  (consisting  of  transla¬ 
tional,  rotational  and  vibrational  energy  of  the  N02  gas)  gets 
effectively  channeled  into  translational  energy  The  observed 
transient  velocity  of  1  63  km  •  s-1  would  then  be  consistent  with  a 
chamber  wall  temperature  of  —  1600  K.  At  this  temperature,  the 
chamber  pressure  would  be  —60  atmospheres  if  the  chamber  had  no 
escape  aperture. 

The  less  volatile  MMH  appears  more  like  a  protracted  flow  of 
liquid  than  a  compact  cloud.  Because  of  the  lower  vapor  pressure,  it 
can  be  expected  that  the  pressure  due  to  evaporation  of  the  fuel  is 
significantly  lower  than  for  the  oxidizer,  which  is  consistent  with  the 
significantly  lower  velocity  measured  for  the  leading  edge  of 
transient  3. 

Evolution  of  Condensates  Injected  into  Space 

In  the  present  section,  the  transient  condensates  are  modeled  in  a 
more  quantitative  way.  Focusing  on  transient  2,  the  decline  in 
integrated  cloud  intensity  as  a  function  of  time  with  respect  to  engine 
shutdown,  shown  in  Fig.  6,  is  indicative  that  the  oxidizercondensates 
are  losing  mass  and  shrinking  in  radius  through  sublimation.  This  is 
analogous  to  the  observations  made  during  space  shuttle  water 
releases  [13].  As  described  by  Kofsky  et  al.  [13],  the  size  and 
temperature  evolution  of  a  droplet  in  space  can  be  obtained  by 
integrating  the  heat  balance  equation: 

47T 

Q(r.  T)  =  y  rVCpCDT  =  Qsll„(r)  +  ()Emh(r) 

~ QUr.T)-Qiub(r,T )  (5) 

where  Q(r,  T)  is  the  particle  heat  gain  or  loss  rate,  which  is  given  by 
the  sum  of  heating  by  solar  and  Earthshine  radiation  and  the  cooling 
through  radiative  and  evaporative  or  sublimation  cooling;  p  = 
1.52  gem-3  is  the  density;  and  Cp{T)  is  the  temperature-dependent 
specific  heat  (constant  pressure)  of  N204.  The  heating  terms  are 
given  by  the  solar  and  Earthshine  absorption  expressions: 

Q(r)  ~  ^(r,  u)Qr2/5(v,  T^dv  (6) 

where  Q  is  the  solid  angle  subtended  by  the  Earth  (— 1.4;rsr  at 
335  km)  or  sun  (6  8  x  10"5  sr),  and  P(v,  Tt)  is  the  Planck  radiation 
function  with  respect  to  the  frequency  v,  where  7,  =  280  and  5800  K 
for  the  Earth  and  sun,  respectively.  The  radiative  loss  term  is  given  by 

Qno(r,T)  =  4nf  e(r,  v)r*P(v,  T,)dv  (7) 

and  the  sublimation  term  is  calculated  from 

Q  subCr)  =  47rrJpAH,obr  (8) 

where  A//Sllb  =  590  J/g  is  the  average  heat  of  sublimation  over  the 
pertinent  temperature  range.  The  emissiviues  were  calculated  from 
the  indices  of  refraction  reported  by  Roux  and  Woodward  [  1 5]  using 
standard  Mie  theory  [16].  Although  Roux  and  Woodward  [15]  do 
not  report  indices  below  wavelengths  of  —2  /an,  we  assume  no 
significant  variation  of  the  real  index  of  refraction  at  shorter 
wavelengths,  and  we  estimate  the  complex  index  of  refraction  from 
the  absorption  cross  section  of  N02  at  600  nm.  This  leads  to  an 
estimated  refractive  index  of  1.53  ±0.005/  at  the  visible  wave¬ 
lengths  of  observation.  Equation  (5)  is  then  coupled  with  the 
expression  for  change  in  radius  by  sublimation  as  expressed  through 
the  temperature-dependent  vapor  pressure  [13].  Parameters  for  both 
N204  and  MMH  used  for  the  calculation  of  r(t )  and  T (f)  arc  shown  in 
Table  4. 

The  bottom  two  charts  in  Fig.  10  show  the  time  evolution  of  the 
diameter  D=2r  and  temperature  of  the  original  1-,  0  1%  and 
0.02-/zm-diam  N204  particles  exposed  to  sunlight  at  an  altitude  of 
—335  km.  The  diameter  is  plotted  as  its  ratio  with  respect  to  the 
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Table  4  Parameters  used  for  fuel  and  oxidizer  particle  evolution  simulations 


n204 

CHjNjHj  (MMH) 

p 

1.52 

0.88 

cp 

-0.11  +  1.23  x  10-27  +  5.58  x  10-5T2  +  1.09  x  10-7T3 

2.74  -  7.24  x  10~ 57"  +  2.48  x  10"672 

A 

590 

875 

log.o/’v.pfn 

9.82-23327-'  +  8  46  x  lO^”2 

7.11  —  1104  67—'  +  1  52  x  10'7— 2 

initial  diameter  D0.  The  initial  temperature  was  chosen  to  be  250  K. 
This  is  near  the  melting  point  of  N204,  which  depends  significantly 
on  the  content  of  NO  Our  calculations  show  that  it  takes  less  than 
1  ms  to  dissipate  the  N204  heat  of  fusion  of  14.7  kJ/mol,  and  so 
neglecting  the  phase  transition  in  the  present  calculations  only 
introduces  a  minor  error.  The  charts  demonstrate  that  the  initial  stage 
is  dominated  by  sublimation,  leading  to  very  rapid  cooling  (within 
~20  ms)  to  below  150  K.  Dunng  this  period,  the  diameter  of  the 
particle  significantly  decreases  below  D/D0  =  0  7.  After  the  first 
~200  ms,  the  temperature  and  size  of  the  particles  declines  only 
slowly  as  the  particle  reaches  a  steady  state  in  which  the  dominant 
cooling  mechanism  is  radiative  and  the  predominant  heating 
mechanism  is  Earthshine  absorption 

The  top  chart  in  Fig.  10  plots  the  predicted  relative  scattered 
sunlight  intensity  as  a  function  of  time  for  the  subject  particle  at  the 
scattering  angle  0=112  deg  (<p  =  180  deg)  that  was  effective  during 
bum  1 .  The  intensities  for  the  different  sizes  are  brought  to  the  same 
scale  by  normalizing  with  respect  to  the  intensity  at  0.3  s,  which  is 
approximately  the  time  after  shutdown  at  which  the  oxidizer  cloud 
intensity  can  be  determined  The  scattered  intensity  from  a  single 
particle  is  calculated  from 

%)=  f  dvg(v)Fs(v)a(v,r,(fi)  (9) 

JO 

where  £(v)  incorporates  the  spectral  atmospheric  transmission  and 
sensor  sensitivity,  Fs(v)  is  the  solar  flux,  and  rr(v,  r,  </?)  is  the  Mie- 
scattering  differential  cross  section.  The  result  is  compared  with  the 
integrated  intensity  evolution  for  bum  1  following  scaling  of  the 
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Fig.  10  Evolution  of  oxidizer  particle  size,  temperature  and  intensity 
and  comparison  with  burn  1  intensity  data. 


experimental  data.  Although  the  observed  decline  in  intensity  with 
time  is  recovered  by  the  calculations,  it  is  seen  that  the  observed 
intensity  declines  more  rapidly  than  the  3  calculated  curves.  The 
calculations  indicate  an  increasing  decline  as  the  particle  size  is 
lowered  from  1  to  0. 1  (i m,  despite  a  decreasing  decline  in  diameter 
This  is  indicative  of  the  transitional  scattering  regime  between 
geometric  and  Rayleigh  scattenng  at  the  respective  particle  sizes. 
The  geometric  scattenng  intensities  are  proportional  to  r2,  whereas  in 
the  Rayleigh  regime  they  are  proportional  to  ~r4.  The  decline  in 
intensity,  however,  is  lowest  for  the  smallest  particle  size  of 
0.02  fj, m  This  can  be  attnbuted  to  the  very  small  decrease  in  size  due 
to  the  very  low  evaporation  rates  at  the  low  near-steady-state 
temperature  of  ~90  K  of  these  particles 

The  significant  pressures  produced  by  the  evaporating  oxidizer 
after  shutdown,  discussed  in  the  previous  section,  raise  the 
possibility  that  the  oxidizer  transient  is  due  to  submicron-sized 
condensates  formed  in  the  expansion  of  the  gas.  It  is  therefore 
possible  that  the  initial  increase  or  near-constant  intensity  seen  in 
Fig.  6  may  be  due  to  continued  condensate  growth  up  to  ^0.3  s  after 
shutdown.  This  explanation  implies  that  the  accommodating  gas- 
condensate  collision  rate  exceeds  the  rate  of  evaporative  loss. 

The  present  calculations  must  be  considered  highly  approximate 
for  several  reasons:  first  and  foremost  is  the  assumption  of  a  single 
particle  size.  Furthermore,  the  lack  of  absorption  coefficients  for 
solid  N204  at  wavelengths  exceeding  2  (im  required  an  estimate, 
which  in  itself  could  explain  the  discrepancy  between  observation 
and  experiment.  Finally,  the  calculations  assume  particles  consisting 
of  smooth  spheres.  As  discussed  by  Kofsky  et  al.  [13],  simple  Mie- 
scattenng  calculations  such  as  the  one  presented  here  do  not  account 
for  particle  albedo  reduction  due  to  surface  roughening  that  may  be 
caused  by  sublimation  This  could  also  account  for  the  more  rapid 
decline  in  intensity  observed  in  the  experiments 

Condensate  Deposition  Rates 

The  STS-1 15  observations  reported  here  have  allowed  the  first 
determination  of  the  angular  distribution  of  condensates  in  the 
exhaust  of  a  PRCS  engine  Next,  we  relate  these  angular  distributions 
to  previous  centerline  deposition  rates  and  suggest  values  for  the 
deposition  rate  as  a  function  of  distance  and  angle.  Previous 
experiments  have  quantified  deposition  rates  on  the  thrust  centerline 
[1,2],  In  the  work  by  Soares  et  al  [2],  a  deposition  rate  of  1.282  x 
10-5  g  *  citT2  •  s-1  was  derived  from  the  study  of  20  80  ms  PRCS 
engine  firings  at  a  distance  of  10.58  m.  The  deposition  rate  was 
denved  assuming  a  total  engine  firing  duration  of  1 .6  s.  From  the 
Brook  equation  for  the  engine-exhaust  number  densities  and  the 
parameters  in  Table  3,  we  can  determine  an  estimate  for  the  total 
centerline  exhaust  particle  flux  F  at  a  distance  of  10.58  m  from 

F  =  N(Re)ve  (10) 

where  ve  =  3.5  km  •  s-1  is  the  far-field  exhaust  velocity,  and  N(Re) 
is  the  Brook  number  density  at  the  axial  distance  Rc  [9,10]. 
Assuming  an  average  molecular  mass  of  23  atomic  mass  units  based 
on  the  known  exhaust  constituent  distribution  [  10],  a  mass  flux  F  •  m 
of  1.56  x  10-3  g  •  cm-2  •  s-1  is  denved.  Thus,  the  mass  flux 
associated  with  the  deposition  registered  by  the  quartz  crystal 
microbalance  (QCM)  in  the  study  by  Soares  et  al  [2]  represents 
approximately  0.8%  of  the  exhaust  mass  flux. 

This  suggests  that  the  QCM  deposition  may  be  associated  with  the 
presently  discussed  transient  condensate  fluxes.  If  this  were  the  case, 
the  deposition  rate  relates  to  the  number  of  engine  bums  instead  of 
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the  bum  lengths  and  amounts  to  approximately  1  fig  •  cm-2  per  bum 
at  the  QCM.  By  relating  the  derived  condensate  mass  flux  fraction  at 
the  QCM  to  the  total  exhaust  mass  ejected  during  a  single  0.08  s  bum, 
which  can  be  derived  from  the  known  PRCS  engine  mass  flow  of 
1 .4  kg  ■  s“\  the  QCM  measurements  suggest  that  1  g  of  condensates 
survive  per  PRCS  bum  at  a  nozzle  distance  of  1 0.58  m.  This  quantity 
is  consistent  with  the  propellant  mass  in  the  engine  vacuum  manifold 
at  shutdown.  Given  the  previously  mentioned  manifold  volumes  of 
—5  cm-3  for  the  oxidizer  and  ^8  cm'3  for  the  fuel  and  given  the 
densities  of  1.52  g  •  cm-3  and  0.88  g  •  cm"3  for  N204  and  MMH, 
respectively,  7.6  and  7  g  of  unpressunzed  oxidizer  and  fuel,  respec¬ 
tively,  are  in  the  vacuum  manifold  at  engine  shutdown  Thus,  if 
the  QCM  signal  is  entirely  due  to  the  observed  transients,  less  than 
10%  of  the  shutdown  transient  material  is  in  condensed  form  at  a 
distance  of  -'-10  m,  assuming  that  the  QCM  measurements 
efficiently  detected  it.  The  condensate  evaporation  studies  applied  to 
oxidizer  particles  discussed  in  the  previous  subsection  indicates  that 
particles  with  radii  between  0.01  and  10  fim  lose  approximately 
60-70%  of  their  mass  due  to  evaporation  after  a  flight  time 
corresponding  to  a  flight  distance  of  1 0.58  m  at  the  average  oxidizer- 
transient  velocity.  Calculations  conducted  for  MMH  fuel,  on  the 
other  hand,  indicate  that  only  approximately  30%  mass  loss  occurs 
during  the  same  flight  distance 

The  present  estimates  based  on  the  QCM  results  must  be  regarded 
as  lower  bounds,  given  the  relatively  high  ejection  velocities, 
particularly  of  the  oxidizer  transient.  The  associated  hypervelocity 
impacts  at  the  QCM  could  result  in  considerable  ejection  of  material, 
resulting  in  lower  readings.  It  can  be  expected  that  the  QCM 
detects  substantially  more  fuel  than  oxidizer  because  of  the 
significantly  lower  velocities  of  the  fuel  transient,  which  we  have 
only  determined  for  the  first  burst  and  which  can  be  expected  to  be 
lower  during  the  persistent  “dribble”  stages  of  the  transient  An  x-ray 
photoelectron  spectroscopy  analysis  of  impact  craters  due  to 
condensates  in  the  PRCS  exhaust  indicated  high-oxidation-state 
nitrogen,  possibly  attributable  to  MMH-nitrate  [2],  thereby 
suggesting  effects  due  to  both  fuel  and  oxidizer.  The  lack  of 
information  on  the  size  distributions  of  oxidizer  and  fuel  particles 
prevents  quantifying  the  optical  signals  in  terms  of  total  condensate 
mass.  However,  the  time  dependence  of  the  GLO  intensities 
associated  with  condensate  solar  scatter  and  the  evolution  of 
transients  observed  in  the  STS- 1 1 5  MAUI  data  indicate  that  the  time- 
integrated  solar- scattered  signal  due  to  fuel  is  more  than  a  factor  of 

3  more  intense  than  that  of  the  oxidizer,  which  appears  only  as  a 
narrow,  compact  feature 

A  Brook  parameter  for  mass  deposition  as  a  function  of  nozzle 
distance  and  angle  is  provided  in  Table  3  based  on  the  centerline 
results  of  Soares  et  al.  [2]  of  1  fig’  cm’2  per  bum  at  10  58  m  and  the 
assumption  that  the  results  reported  by  Soares  et  al.  are  directly 
related  to  the  observed  transients  2  and  3.  The  Brook  model  does  not 
take  the  high-angle  tail  of  transient  3  into  account  We  therefore  also 
provide  suggested  parameters  for  the  Lorentz-type  distribution, 
assuming  that  the  condensate  deposition  is  primarily  due  to  the  fuel 
(transient  3).  The  high-angle  tail  may  be  a  consequence  of  fuel  film 
cooling  of  the  nozzle,  causing  liquid  ejection  at  larger  angles. 

Soares  et  al  [2]  observed  and  quantified  craters  formed  on 
aluminum  and  Kapton  samples  exposed  to  84  PRCS  and  17  vernier 
reaction-control-system  engine  firings.  Most  of  the  crater  diameters 
on  aluminum  were  less  than  10  /imin  diameter.  Studies  that  relate 
crater  size  to  particle  diameter  find  that  at  impact  velocities  of  around 
2  km  ■  s"1,  craters  tend  to  be  similar  to  or  less  than  the  particle 
diameter,  whereas  at  1  km  •  s"\  the  crater  diameters  drop  signifi¬ 
cantly  below  the  particle  diameter  [17,18].  In  the  case  of  the 
aluminum  sample,  the  most  abundant  crater  sizes  were  less  than 

4  fim.  Relating  the  number  of  such  craters,  449  mm-2,  to  similarly 
sized  particles  leads  to  the  conclusion  that  the  mass  deposition 
associated  with  the  craters  is  approximately  2  orders  of  magnitude 
smaller  than  that  detected  by  the  QCM  without  taking  material 
ejection  into  account.  It  is  possible  that  the  fast-exhaust  species 
associated  with  transient  1  may  be  primarily  responsible  for  the 
cratering.  Unfortunately,  that  transient  was  too  weak  and  its 
observation  time  was  too  short  to  allow  a  meaningful  analysis. 


Conclusions 

Imagery  of  PRCS  engine  far-field  exhaust  was  analyzed  for  22 
attitude-control  firings  during  an  STS- 1 1 5  terminator  pass  over  the 
Maui  Space  Surveillance  Site  (MSSS).  This  imagery  was  interpreted 
in  terms  of  sunlight  scattering  of  exhaust  condensates.  Each  firing  is 
associated  with  three  transients:  one  weak  transient  at  engine  startup 
and  two  intense  transients  at  engine  shutdown.  An  analysis  of  the 
transient  velocities  finds  that  the  startup  transient  is  close  to  the 
freestream  exhaust  velocity  of  3  km  -  s-1,  and  the  two  shutdown 
transients  have  average  velocities  of  1.63  and  0.7 1  km  •  s“l,  respec¬ 
tively.  Although  the  startup  transient  is  attributed  to  initially 
unbumed  propellant  in  either  a  condensed  or  gaseous  state,  the 
shutdown  transients  are  attributed  to  oxidizer  and  fuel  solid  particles. 
The  compact  shutdown  oxidizer  transient  appears  before  the  fuel 
transient  due  to  its  higher  vapor  pressure.  The  shutdown  fuel 
transient,  on  the  other  hand,  exhibits  an  initial  burst  followed  by  a 
persistent  dribble  trail.  The  time  dependence  of  the  shutdown 
oxidizer  transient’s  integrated  intensity  suggests  that  particles  of 
approximately  1  fim  diameter  or  smaller  are  the  source  of  the 
observed  solar  scattering.  From  the  imagery,  the  angular  distribution 
of  the  PRCS  exhaust  condensates  is  derived.  An  analysis  of  the 
angular  distribution  of  the  shutdown  oxidizer  transient  demonstrates 
that  it  is  well  described  by  a  Brook  source  distribution  with  an 
angular  spread  slightly  narrower  than  that  specified  for  the  gaseous 
exhaust.  A  similar  analysis  on  the  initial  burst  of  the  shutdown  fuel 
transient  indicates  that  a  Brook  distribution  does  not  capture 
considerable  material  at  large  angles. 

The  present  analysis  in  context  with  axial  deposition  measure¬ 
ments  [2]  of  a  shuttle  PRCS  engine  in  space  suggests  that  at  least  1 0% 
of  the  propellant  manifold  volume  is  ejected  as  condensates  at 
shutdown.  This  assumes  that  the  deposition  signals  are  caused  by 
the  shutdown  transients.  The  higher  vapor  pressure  and  associated 
sublimation  rates  of  the  oxidizer  suggest  that  condensate  release  and 
deposition  are  more  important  for  the  MMH  fuel.  This  is  consistent 
with  spacebome  measurements  using  the  GLO  instrument,  for  which 
the  effects  due  to  condensate  solar  scatter  are  observed  beyond  10  s 
following  a  PRCS  bum.  The  extended  solar  scattering  signals  are 
attributed  to  the  persistent  MMH  dribble. 
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